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The Phase Diagram of QCD 
M.A. Stephanov'^ 
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We show that current experimental knowledge of QCD together with general model 
independent arguments such as continuity, universality and thermodynamic relations, as 
well as the information gained from various models can be used to constrain the phase 
diagram of QCD as a function of temperature and baryon chemical potential. 



1. Introduction 
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Q I Understanding the phase diagram of QCD as a function of temperature, T, and chemical 
potential of the baryon charge, /i, is an important ingredient in the analysis of the results of 
heavy ion collision experiments. We know a good deal about the behavior of QCD at finite 
temperature. The basis for our understanding is provided by the first principle lattice 
QCD Monte Carlo calculations. In particular, we know that QCD with two massless 



D i quarks undergoes a phase transition to the quark-gluon plasma phase at a temperature of 
I about 160 MeV. In contrast to that our understanding of the behavior of QCD at finite 
baryon charge density, or finite /i, is extremely poor. This is due to the unfortunate fact 
^ ■ that lattice QCD calculations based on Monte Carlo methods are not possible because 
the measure of the Euclidean path integral is complex when /i is not zero. In this report, 
based on the work we put together available experimental facts about QCD, results 
from various models, and apply general model independent arguments such as continuity, 
universality and thermodynamic relations in order to construct the phase diagram of QCD 
in the T/i plane. 

We perform our analysis for the two-flavor QCD — a well-known and phenomeno- 
logically successful approximation to real QCD. The effects of the strange quark and 
electroweak interactions lead to quantitative and in some cases qualitative modifications 
which we also consider. The theory is described by a partition function: 

Z = e-"(^'^)/^ = J VAV^jViPexpi-Ss} ■ (1) 

The Euclidean action, Se, is given by 



Se = J dxo I 



1 ^/ / 



(2) 



where Nf = 2 is the number of flavors, Nc = 3 is the number of colors, and mj = m = 
is the quark mass. The Euclidean matrices 7^ are hermitian. The normalization of /j, 
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differs from tfie normalization customary in lattice calculations by a factor l/N^ (i.e., the 
baryon charge of a quark). 

What do we actually mean by understanding the phase diagram? The most prominent 
features of a phase diagram are phase transitions. They manifest themselves through the 
singularities or discontinuities in the dependence of various thermodynamic observables 
on the parameters T and /x. Such observables can be obtained by differentiating the 
thermodynamic potential f2 with respect to T and /i: 

- dVt dVt - dVt 

nV = j:{i:noi^f) = - — - sV = - — - and also mNfV = -g^; (3) 

where n is the baryon number density, and s is the entropy density. These are the densities 
of extensive quantities, such as baryon charge and entropy, per volume, V. The pressure, 
pV = —dQ/dV = —Q, is not independent from T and fi (and also m): 

dp = sdT + nd^ + {il)ip)N fdm. (4) 

This equation can be used to derive Clapeyron-Clausius-type relations between the slopes 
of the first-order transition lines and the discontinuities of s, n and {ipip) 

There are two thermodynamic observables which turn out to be more useful than oth- 
ers in discovering phase transitions: n and This is because both are good order 
parameters, i.e., they vanish identically in one phase and are nonzero in the other. There 
must be a singularity, and thus a phase transition, separating such two phases. Theoret- 
ically, we understand, at least qualitatively, the behavior of {^ip) because it is an order 
parameter of a global symmetry, SU(2)iXSU(2)ji;, and it distinguishes two phases with 
two different realizations of this symmetry: spontaneously broken and exact. However, 
phenomenologically we know little about this phase boundary (yet!). On the contrary, it 
is harder to understand theoretically the behavior of n, but, fortunately, we happen to 
live near the phase boundary separating phases with = and n 7^ 0. As a result, we 
have a good empirical and quantitative knowledge about this phase transition. 

2. Zero T 

To see why n can serve as a good order parameter consider the partition function (|l|) 
written in the form of the Gibbs sum over the quantum states of the system characterized 
by their energy, and baryon charge, N: 

In the hmit T — > 0, the state with the lowest value of — fiNa determines the properties 
of the system. This is the ground state at given /i. Let us introduce 

Ho = mm{Ea/Na) ■ (6) 

As long as /i < yUo no state with nonzero can compete with the vacuum state [E = 0, 
A^ = 0) for the role of the ground state. Therefore, as long as T = and /i < fiQ, the 
equation n{fj,) = holds exactly. If n 7^ when > /xq, the point /i = /io must be a 
singular point. 
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Figure 1. Schematic dependence of the baryon charge density on the chemical potential 
at T = (a) in QCD (/xq ^ mjv - 16 MeV) and (b) in QCD+ (/^o ^ thn - 8 MeV). At 
III the chiral symmetry restoration transition occurs. 



In QCD without electromagnetism this singularity is a first order phase transition 
separating vacuum phase from nuclear matter phase, distinguished (at T = only) by 
the order parameter n. The function n(/i) has a step at /i = /io ~ — 16 MeV equal to 
the density of nuclear matter no ~ 0.16 fm~^ at zero pressure as in Fig. |l|a. The slope of 
n(/i) just above /io is also known. In the real world, the electromagnetic interaction and 
the presence of electrons produce a tiny step at /i = /xo~mAr — 8 MeV to the density of 
iron. The step to the density of neutron matter occurs at somewhat larger value of /i (see 
Fig. |l]b) at non-zero pressure (as in neutron stars). 

Increasing fi further takes QCD into the region about which we have very little reliable 
theoretical or experimental information. Various interesting phenomena have been pre- 
dicted using different models. Rigorously, however, one can only expect that the function 
n(/x) continues to grow, to satisfy thermodynamic stability. Since n is non-zero it cannot 
serve as a good order parameter anymore and we shall turn to another one: (ipip) ■ 

At very large /i the ground state of the system is to a good approximation a Fermi 
sea of quarks. (This fact also leads to an interesting possibility of quark-quark pairing 
around the Fermi surface and color superconductivity 0.) Thanks to the asymptotic 
freedom and the screening of color interactions by the sea of quarks, nonperturbative 
effects are suppressed. This motivates the conclusion that at very large fi the condensate 
(ipip) vanishes (provided that quark masses are zero). If we denote by fii the value of 
II such that {ipip) = for /i > /ii and {ipip) ^ otherwise, then the point n = Hi 
must be a singular point. It separates two phases with distinct realizations of the global 
SU(2)LxSU(2)i{ chiral symmetry. 

There has been a wealth of theoretical research on the phase transition at /x = /ii (see, 
e.g. 1^). However, since no first principle lattice calculations are possible at present, the 
best we have is a collection of estimates obtained in various models approximating the 
behavior of QCD at the chiral symmetry restoration transition. The common denominator 
seems to be that the transition is most likely of the first order. The value of iii is 
somewhere of the order of 1 GeV. Also, the empirical fact that stable nuclear matter with 
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n 7^ and ((ipip) 7^ exists indicates that /ii > /iq. (The strange quark tends to lower the 
value of /ii- An interesting, but empirically disfavored, possibility that /ii < /io, strange 
quark matter at zero pressure, arises in this case 0].) 

3. Nonzero T 

At finite T the baryon density n is no longer a good order parameter. However, since 
the transition at /iq is of the first order, continuity ensures that it remains first-order 
for sufficiently small T. The end-point of this transition is a critical point in the Ising 
universality class, which is probed in the multifragmentation experiments 0. 

The chiral condensate {ipip) is a good order parameter at T 7^ 0, as long as m = 0. This 
means that the regions {ipil') 7^ at low T and and {^Jip) = at high T and ^ must 
be separated by a phase transition. At yU = lattice simulations predict that the chiral 
symmetry is restored at Tc ~ 160 MeV for two- flavor QCD ||^. A beautiful argument 
suggests that this transition is likely to be a second order transition in the universality 
class of 0(4) spin models in 3 dimensions. 

At nonzero /i there is a line of second order phase transitions in the 0(4) universal- 
ity class (the /i direction is not relevant near this critical point). Since this line cannot 
terminate, and since at T = the transition is, presumably, of the first order, a log- 
ical possibility arises that the transition turns first order in some point T3, /^s on the 
phase diagram. This point has been observed in various models of QCD chiral phase 
transition |^^,|T|,|TD| . We wish to point out that the critical behavior near this point is 



determined by universality (the observation also made independently in |jTO[)- Study of 
tricritical points shows that the upper critical dimensionality for such a point is 3 
and, therefore, in QCD with two massless quarks the critical behavior near the tricritical 
point must be given by the mean field theory up to logarithmic corrections. 

In this work we used a random matrix model to describe the chiral phase transition. 
In accordance with generic expectations and other models it predicts a tricritical point. 

4. A random matrix model at finite T and jjL 

A successful and very simple model which describes the degrees of freedom of QCD 
related to the spontaneous breaking of chiral symmetry is the random matrix model. It 
is based on the famous observation of Banks and Casher that the value of {^ip) is related 
to the density of small eigenvalues of the Dirac operator pev(O). Since we do not need to 
describe all the degrees of freedom of QCD, but only those relevant to chiral symmetry 
breaking, a natural and simple approximation of the Dirac operator by a random matrix 
arises |T^ . This approach is reminiscent of the one introduced by Wigner in the study of 
spectra of heavy nuclei. The power of the random matrix model in describing the chiral 
symmetry breaking and restoration at finite chemical potential is in the fact that, on the 
one hand, this model shares an important property of QCD at /i 7^ — the complex 
fermion determinant, and it is exactly solvable, on the other hand. A successful example 
is the explicit demonstration that at nonzero quenched QCD is not a smooth Nf 



limit of real QCD |T3|. 



The phase diagram calculated in the random matrix model IQ is shown in Fig. ^. We 
see that the second order phase transition line in the plane m = turns into a first order 
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Figure 2. Phase diagram of QCD with two hght flavors of mass m as calculated from the 
random matrix model. The almost parallel curves on the wing surface are cross sections 
of this surface with m =const planes. The units of m are 100 MeV, of T are Tc ~ 160 
MeV, of /i are /ii/0.53 ~ 2300 MeV, with the choices of Tc and /ii from the text. 



line at the tricritical point. The coordinates of this point are given by: T^/T^ ~ 0.78, 
li^/Hi « 0.61. Taking = 160 MeV and /ii = 1200 MeV, we find that T3 ^ 120 MeV and 



/i3 700 MeV. A similar estimate has been obtained recently using a different model [10 



5. Conclusions 

We conclude with a sketch of the phase diagram (Fig. ^ of QCD with two massless 
quarks which we find by analyzing the behavior of two thermodynamic quantities: n and 
(ipip). These quantities are distinguished by the fact that both are good order parameters 
in a certain sense: they identically vanish in some region of the phase diagram and are 
non-zero in the other. The change of the behavior of such a parameter from one region 
to the other is qualitative and must proceed through a thermodynamic singularity. 

Perhaps, the most interesting feature of this phase diagram is the presence of the 
tricritical point. At small nonzero quark mass the main change in the phase diagram is 
the disappearance of the second order phase transition line (see Fig. The first-order 
line remains, but it no longer separates phases with different symmetry properties — the 
chiral symmetry is explicitly broken. What is important is that the criticality at the 
end-point of the first-order phase transition line remains. 

The strategy for locating this end-point in the heavy ion collision experiments is dis- 



cussed in [14|. The signatures proposed in [14] are based on universal thermodynamic 



properties of the critical point. For example, the divergence of the heat capacity will lead 
to suppression of the event-by-event fluctuations of the apparent temperature. Another 
signature is due to the long-wavelength fluctuations of the sigma field near the critical 
point, which leads to enhanced production of soft pions. An inspiring example that the 
study of critical behavior in heavy ion collisions may, in principle, be possible is provided 
by the multifragmentation experiments [0, which can probe the end-point of the nuclear 
liquid-gas phase transition. 
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Figure 3. A schematic phase diagram of QCD with 2 massless quark flavors. Other phase 
transition hues are possible, for example, in the low temperature region to the right of 
yUo- Another example is a transition associated with color superconductivity plotted as a 
dashed line. Thicker lines are first-order phase transitions. The — T3 line is a second- 
order phase transition. The tricritical point is at T3, fi^ and the critical point of the 
nuclear matter liquid- gas transition is at Tq. 
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